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Abstract: We consider an optimal control problem for the three-dimensional 
non-linear Primitive Equations of the ocean in a vertically bounded and hor- 
izontally periodic domain. The observation operator maps a solution of the 
Primitive Equations to the trajectory of a Lagrangian particle. This paper 
proves the existence of an optimal control for the regularized problem. To do 
that, we prove also new energy estimates for the Primitive Equations, thanks 
to well-chosen functional spaces, which distinguish the vertical dimension from 
the horizontal ones. We illustrate the result with a numerical experiment. 
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Optimal control of the Primitive Equations of 
the Ocean with Lagrangian observations 

Resume : On considcrc un probleme de controlc optimal pour Ics equations 
primitives non lineaires de I'ocean en dimension trois d'espace, dans un domaine 
verticalement borne et horizontalement periodique. L'operateur d'observation 
associe a une solution des equations primitives la trajectoire d'une particule 
lagrangienne transportee par le flot a profondeur donnee. Ce travail montre 
I'existence d'un controle optimal, apres regularisation du probleme. Pour cela, 
on etablit en particulier de nouvelles estimations d'energie pour les equations pri- 
mitives dans des espaces fonctionnels bien choisis, qui distinguent la dimension 
verticale des dimensions horizontales. On illustre le resultat par une experience 
numerique. 

Mots-cles : controle optimal, equations aux derivees partielles, equations 
primitives, simulation numerique 
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Introduction 

The ocean plays a major role in governing the earth climate. Physical oceanogra- 
phers and climatologists work toward a better knowledge of the ocean properties 
(currents, temperature, salinity, marine biology, etc.). Some mathematical tools 
are involved in this progress, in particular Data Assimilation methods. Data 
Assimilation covers all mathematical methods which allow to blend optimally 
all sources of information about the ocean (measures, model equations, errors 
statistics) in order to improve ocean modeling, forecasts or climatology. One 
class of data assimilation methods, called Variational Data Assimilation [HIS], 
is based on optimal control theory [B]. The idea is to compute numerically the 
solution of an optimal control problem, in which the cost function represents 
the misfit between the observations and their model counterpart. 
This paper deals with a particular problem of data assimilation for the ocean, 
namely the assimilation of Lagrangian data. The ocean is mainly observed at 
the surface, thanks to observing satellites. In-situ data are sparsely sampled in 
time and space, and it is therefore important to make the most of their infor- 
mation. Lagrangian data consist of positions of floats drifting at depth (around 
one thousand meters deep), they give information about in-depth currents. The 
problem of variational assimilation of Lagrangian data has been studied numer- 
ically in [TD]. 

In this paper we investigate the theoretical justification of this problem. We 
prove the existence of an optimal control for Lagrangian observations for the 
Primitive Equations (PEs) of the ocean. To do so, we had to establish local 
existence and unicity of strong solutions for the PEs, allowing existence of La- 
grangian trajectories, so that the problem can be formulated. The problem 
of local existence and unicity of strong solution for the PEs has been studied 
by Lions, Temam, Wang and Ziane |H1 [T3]. Due to the dissymmetry between 
the vertical dimension and the horizontal ones in our chosen domain (vertically 
bounded and horizontally periodic) and also in the PEs (as the equation for the 
vertical velocity is degenerated, contrary to Navier-Stokes equations), we had 
to introduce new functional spaces. In these spaces we successively prove new 
energy estimates for the linear PEs and existence of strong solutions for the 
non-linear PEs. 

This paper is organized as follows. In section [l] we state the equations, the 
functional spaces, the cost function, the optimal control problem and the main 
results of the article. In section [2] we prove new energy estimates for the linear 
PEs and local existence and unicity of the non- linear PEs. In section [3] we prove 
the existence of an optimal control. Finally, in section [4] we present a numerical 
illustration of this problem. 
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1 Statement of the problem and main results 

1.1 The Primitive Equations of the ocean 

We consider the Primitive Equations of the ocean in a three dimensional domain 
(see UM)- 

dtu — vAu + {U.V2)u + wd^u — av + d^p = in fi x (0, T) 

dfV — vAv + (C/.V2)u + wdzV + au + dyp — 

d^p -(36 = 

dt9-iyA0+{U.\72)O + wd,e + jw = O mnx{0,T) (1) 

w{x, y,z) = - J^ dj;u{x, y, z') + dyv{x, y, z') dz' in Vl x (0, T) 

U{t = 0) = t/o, 0{t = 0) = 6*0 in f} 

with 

- f^ = T2 X (0, a) the physical domain, with T^ = (M/27rZ)2 the bidimcnsional 
torus, such that the domain is periodic in the horizontal directions x and y, 
vertically bounded in z, with fixed depth; 

- (0, T) the time interval; 

- U = {u, v) the horizontal velocity vector, w the vertical velocity, 9 the 
temperature and p the pressure; 

- Uq = {uq,vq) and 9o the initial conditions; 

- V2 = {dx,dy) the horizontal 2D gradient operator, (V2-) the horizontal 
divergence operator, V = {dx,dy,dz) the 3D gradient operator, A = dxx + 
dyy + dzz the 3D Laplacian; 

- a, p, 7, P physical constants. 

The boundary conditions are the following: 

M, V, 9 are periodic in x, y 

u = 0,v = 0,9 = 0onT'^ x{z = 0,z = a}x {0, T) (2) 

XLo ^^" + ^y^ dz = on T2 X (0, T) 

We denote by X{t) = {u{t),v{t),9{t)) the state vector of our system, and by 
Xq = {uo,vo,9o) the initial state, which will be our control. 
We will focus on smooth solutions of the Primitives Equations (PEs), so that 
the cost function (involving Lagrangian trajectories) can be defined. To this 
end we first define the following functional spaces, for rn, G N: 

LlH'Jl = {ue L^{n) periodic in x,y, d^^yU e L^{n),\fa e N'\ \a\ < m} 

U'"'+^ = {X = {u,v,9) e {LlH^yf, periodic in x,y, 
X = 0onT2 X {z = 0,z = a}, 

X2L0 "^^^ ^ '^v^ dz = on T^ 
VX=(Vu,Vz;,V0)e((L2i/™)3)3} 

W"+^ = {w e L^iJ™, periodic in X, J/, 

M = on T^ X {z = 0, z = a}, 
Vu e {LlH^yf } 
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associated to the following scalar product and norms: 
(mi,U2)l2h™= ^ d^yUid^yU2dxdydz 

I a I < m 

{Vui,Vu2)(LiHsi)^ ^ Yl d^yVui.d"yVu2dxdydz 



\oc\<.7n 
{Xi,X2)u^ + i == {ui,U2)lIH^^ + (Vui, V'U2)(L2^™)3 

Hvi.V2)lih^^ +(Vwi,Vi;2)(l2^™)3 

{K is a "large" constant which will be set later) . 

We define on (7i™+^)'^ the same scalar product and norm as on U"^'^^ . 

Remark 1 The functional spaces Zi'"+^ and (7Y™+^)^ are not interpolation 
spaces. In the sequel, m is a fixed integer larger or equal to 2. 

In this framework, we have the 

Theorem 1 Let m > 2 be an integer and X^ = (uo,uo,6'o) G W^^^ . If K is 
large enough, there exists t* > Q with t* — t*(a,/3,7, i^, llXgH^^m+i) and there 
exists a unique solution X{t) — {u{t),v{t),9{t)) of the PEs Hy with boundary 
conditions Q) such that 

XeC([0,r];Z^™+i), dtXeL^{0,t*;LlH^y) 

Moreover, we have: 

2 , -^ / lla Vf„\u2 -,„ ^ ^'-' II V l|2 



X{t)\\i,^^, + - \\dtX{s)\\i„,ds < -\\Xo\\t,„.+, (3) 



I'./o '"" S 



for all t G [0, i*], where 5 depends on t* . 
This result is proven in section [2J 

1.2 The Lagrangian observations and the cost function 

We can assume, without loss of generality, that there is only one drifting fioat, 
and that its position is observed at only one given time ii. Its position ^(i) — 
(^^(t),^^(i)) in the plane z = zq is solution of the following differential equation: 

m = 6 

The following proposition is an easy consequence of theorem IT] 

Proposition 1 Under the hypothesis of proposition [71 the unique solution X 
of the PEs (Up and ^ is continuous in time and in z, Lipschitz in {x,y). 
Moreover, for all ^o G "^^ o.'^'d zq G [0,a], there exists a unique Lagrangian 
trajectory, solution of equation (Op, associated to X, ^q and zq. 
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We then define the following cost function: 

j{x,) = iu{t,)-dr + iiixoi|2,,„+, 

J°{Xq) + ujJ\Xq) 

with: 

- d= ((^1,(^2) the observation; 

- 771 an integer, to > 2, w a positive constant; 

- 1 1 . 1 1 the Euclidian norm in the 2D plane z — z^. 
The observation operator is thus defined as follows: 

g(ii;Xo)-C(ti) (6) 

where ^ is defined by equation (H) where the velocity field U = (m, v) is solution 
of the PEs ([TJ and ^ initialized with Xq. 

Remark 2 Contrary to the classical theory 0/ J.-L. Lions ^1, the observation 
operator is non linear; moreover it is defined as a function of either the initial 
state Xq, or as a function of the complete velocity field {U{t),t G [0,ii]}, and 
not only U(ti). 

1.3 Statement of the problem and main result 

The optimal control problem associated to the observation of Lagrangian data 
is the following: 

Problem 1 Let d £ R- be an observation. We look for an optimal control 
Xq G W"^^ solution of the following minimization problem: 

J{X*)= -^nf J{Xo) 



where the cost function is defined by [Ey, the state equation by { 1\2 ) and the 
observations by Q). 

The main result of this paper is the 

Theorem 2 There exists an optimal control Xq E W^^^ solution of problemyn 

This result is proven in section [3] 

2 Existence of strong solutions for the Primitive 
Equations of the ocean 

In this section we prove theorem[T]in three steps: first we prove energy estimates 
for the linear PEs, then we prove estimates for the non linear terms of (111 and 
then we prove theorem [l] 
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2.1 Energy estimates for the linear Primitive Equations 

We consider the following linear Primitive Equations: 

in Q. X (0, T) 



dtu — vlS.u — av + dxP = Fi 
dtv — vAv + au + dyp = F2 

dtO - vM + 7u; = F3 



in r^ X (0, T) 

w{x, y,z) = ~ J^ dxu{x, y, z') + dyv{x, y, z') dz' in fi x (0, T) 

U{t = 0) = C/o, 6{t = 0) = 6*0 in n 

with the same notations as in section [O] and boundary conditions (pi. 
In the sequel we will use the following notations: 



(7) 



/// 

11/11 

||/||2,m 

ll/llm 

ll(/l,/2)||2,,n 

||(/l,/2, /3)||2,m 

The following proposition holds true: 



/o Jq f{t,X;y^ z) dx dy dz dt 

ll/IU2(0) 

WfU- 

ll(/l,/2)||(L2ff™)2 
||(/l,/2,/3)||(L2//™)3 



Proposition 2 For all K large enough, for all T > 0, there exist constants 
Ci{a,i',K,j,l3), C2{K,v), Cz{a,v) and Ci{u) such that, for all X^ e U"^^^ , 
F e L^{0,T; L^H"y), the unique solution X{t) of the linear PEs tw satisfies: 

X{t)eC{[0,T],U"'+^) 



Moreover, the following inequality holds true: 

\\xmi^ + \\vx{t)\\i^ + 1 /J ||a,x(.)|ii,„ ds 

< e^^'{C2\\Xo\\l„, + CsllVXolli,™ + C4 /; ||F(.)||i„ ds) 
for all te [0,T]. 



(8) 



Proof. Classical variational methods (see for example [8j [13]) prove that 
for Xo given in U"^+^ and F e L^{Q,T-LlHj;}y), there exists X{t) at least in 
L2(0,T;V) nC([0,T];H), where V and U are classical spaces (see [H [Tl|3]) 
defined as follows: 

Definition 1 Let 

El = {U = {u,v) eC°°{^Y,u,v periodic in x,y, 
u = 0, -y = on T^ X {z = 0, z = a} 
/„" dxu{x, y, z') + dyv{x, y, z') dz' == 0, V(a;, y) G T^} 

-E2 = {9 E C°"{il.),9 periodic in x,y, 
9 = onT'^ X {z = 0,z = a}} 

Then Tii (respectively Ti-i) is defined to he the closure of Ei in i^(51)^ (resp. 
L'^{fl)), and Vi (resp. V2) is the closure of Ei (resp. E2) in H^{fl)^ (resp. 
H^{n)), and finally H = Hi x H2, V = Vi x V2. 
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Thus it suffices to prove (l8|. To this end, we successively state four energy 
estimates: first an estimate of ||X(t)||i2 and similarly of ||X(i)||i2^m , then 
an estimate of ||VX(i)||i2 and similarly of ||VX(i)||i2//m . 
To obtain energy estimates for ||jr(i)||i2 ^, we multiply equations ([7| by u, v, 
w, K9 and we integrate in space an time. Thus we have: 



(9) 



Ti + T2 + Ta + r4 + T5 = Te 

Ti = // dtuu + dtvv + KdiOe T2 = Jf -vAuu - vAvv - vKMO 
T3 = JJ —avu + auv T4 = J J —I39w + Kjw9 

T5 = // d,pu + dypv + d,pw n = JJ Fiu + F2V + KF39 

We integrate by parts using conditions (l2| and we get: 

Ti = UMm'+Mm'+K\mt)r-\\uor-\\vor~\\eor) 

T2 = vjl\\Vu{s)f + \\Vv{sW + K\\Ve{sWds 

= ^/„*||VX(.)|p (10) 

T3 = 

Ts = - ^^ p{dxU + dyV + dzw) + ^^ j,^2P{w\z=i-w\z=a)dxdydt 
= 



Then ^ and ^ give: 

\\X{t)f+2v f \\\7X{s)f = \\Xo\\^+2 f f {Fiu+F2V+KF39)-2{Kj-(3) f f w9 

(11) 
Then we give a bound of the right hand side of equality (111. First we establish 

the following useful inequality for 1 1 w 1 1 : 

||u;|P = J^\ Jg dxU + dyvdz'l"^ dxdydz 

< J^2zJ^\da:u\^ + \dyv\^dz'dxdydz (12) 

< a^{\\dM\^ + \\dyv\\^) 

Thanks to ([121 we have: 



2jJ{F,u + F2V + KF39) < /„* ||Fi(,s)||2 + ||F2(,s)||2 + /f ||F3(s)|p 

+ \\u{s)\\^ + \\v{s)f + K\\9{s)\\^ ds 

= j;;\\F{s)r + \\x{s)rds 

-2{Kj-(3)JJw9 < \K^-P\J^{\\wr + \\9r) 

< 1X7 - P\ /o a^d^l^ + a^dyvr + \\0r ds 



(13) 
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Then (jUj) and ([13| lead to: 

\\x{t)r+2^j^\\vxis)r < \\xof+j^\\F{s)r+j^\\x{s)rds 

+ J^\\d.ur+\\dM\'ds 

+ \Kj ~ f3\ /o* a^d^u^ + a^dyvf + \\Ofds 

< \\xor+jt\\F{.s)r+j^\\uis)r 

+{K + \Kj - p\)\\0{s)f ds 

+ /„* niax(l, a2|i^7 - /3|)|1 VC/(s)||2 ds 

We proceed similarly to obtain a bound on ||X(t)||i2jLfm . First we note that if 
u, V, w, 9 and p satisfy equation (m), then for all a G N^, i9"yM, d'^yV, d^yW, d"y9 
and d"yP satisfy the same equation (where Fi is replaced by d"yFi and Xq by 
d^Xo) and the same boundary conditions. Therefore we can apply the same 
calculations and obtain the following inequality: 



ll^(i)llL. + 2^/ol|VX(s)|||,„ < ll^ollL„ + /o 11^(^)111™ 

+ Io\\Uis)\\h^ + {K+\K^-P\)\\9is)\\l^ds 
+ />ax(l,a2|ir7-/3|)||Vf/(s)||2„ds 

(14) 
To establish the second type of estimates, we then multiply the equation by dtu, 
dtv, dtw and Kdt9 and we integrate over ft x (0,t): 

Ti + T2 + T3 + r4 + ^5 - Te 

Ti = // dtudtu + dtvdtv + Kdt9dt9 

T2 = // -vAudtu - vAvdtv - vKMdt9 

T3 = JJ —avdtu + audfV 

Ti = // -I39dtw + K-/wdt9 

T5 = // dxpdtu + dypdtv + dzpdtw 

Te = JJ Fidtu + F2dtv + KFzdt9 

Using integration by parts and boundary conditions (pi) we get: 

Ti = Jl\\dtu{sW + \\dtv{s)f + K\\dtu{s)rds 

= JlwdtXisWds 

T2 - f||VX(t)||2-|||VXo|p 

T4 = JJiKj + P)wdt9 + J^^{9oWQ-9{t)w{t))dxdydz 
T5 - 
Thus: 
2j^\\dtX{s)\\^ds + i^\\WX{t)\\^ = i^WWXof + 2 Jjavdtu-audtv 

-2 JJiKj + [3)wdt9 + 2 J^{9{t)w{t) - 9QWa) dx dy dz 
+2 J J Fidtu + F2dtv + KF3dt9 

In order to get a bound on the right hand side we use the following inequality: 

xy<7:x^ + 7^y^ Vx, ?/ e M, Ve > 
2 2e 
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Thus we obtain, for all positive real numbers £i, 2 < i < 5: 

2 If avdtu - audtv < a J^{e2\\dtU{s)\\^ + j^\\U{s)\\^) ds 
-2jJ{Kj + l3)wdte < {Kj + f3) J^ e3\\dt0{s)f ds+ 

2j^{9{t)w{t)-9oWo)dxdydz < a^ei{\\dMtW + Wdy^itW) + rM^W 

2jjF,dtu + F2dtV + KF3dte < J^{±\\F{s)f + e5\\dtX{s)\\^)ds 
And finally we get: 

+ J^{2~ae,-s,)\\dMm'~imt)r 
+{2K - {Kj + /3)e3 - Ke5)\\dte{s)f ds < {i^ + a^)\\WUo\\^ + i^J/||V0o|P + 11^011' 

+ iK^ + (3)Ioi\\yUis)rds + aj;^^J\Uis)rds 
+ JoJ;\\F{s)rds 

As previously we obtain the same result for the derivative in x and y: 



{iy-a^e,)\\VUmhn + K,^\\V9{t)\\ 



2,m 



-ill^Wlli,™ + /o(2 - «£2 - e5)||9tf/(s)|lL„ 
+ {2K~{Kj + (3)e3~Ke5)\\dte{s)\\l^^ds < {j, + a^)\\WUo\\l 

ft „2 



^'^l|veollL„ + ll^o||i, 






The £i are chosen as follows: 



i i 

£2=-, £3=-, 
a 7 


£4 = 


V 

2a2' 


£5 = 


i 

'~ 2 


and thus we get: 










fl|vf/(i)||L„ + ifi^llv0(i)iiL„ 










-'-^\m)\\lra + SlmU{s)\\l„. 











+(f - 7/3) II 9*0(5) II 1,„ ds < {1. + a')\\VUo\\i„. + i^HlV^olll™ + ll^olll. 



+ j;a^j{Kj + P)\\WU{s)\\l„,ds 
+ IoO^'\\Ui-s)\\l^ds + j;;2\\Fis)\\ 



2,m "^ 



We then add equation (14 1 and equation (151 multiplied by - to obtain 



(15) 



iic^(^)iiL„ + iK- ^mmhr, + iivc/(i)ii|„ + 2K\\v0m\ 



+ /o^liatt/(.)||t„ + (f-7/3)^||5,0(.s)||l„d5 + /j2H|Vf/(.)||i,„ + 2^X||V0(.)|| 



2,m 



ds 



< WUoWLn + iK + ^)||eollL„ + (2 + ^)||VC/o||L„ + 2K\\\/eo\\l^ 

+ /„*(l + i)||F(.)||i,„ds/„*(l + 2^)||t/(.)||l,„ds + /„*(i^+|if7-/3|)l|e(s)||i,™ds 
RR n"6490 + /o*(max(l, a2|i^7 - /3|) + ^^(^7 + /3))|| VC/(s)||i,„ ds 
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For K > 2max(^, 27/3), for all T and for all t e [0, T] we have: 

2 illV7\^/'+M|2 I / IIQv/'„M|2 J„ ^ r< \ r< I W V I „\\\'^ illV7V/'„M|2 



^ Jo Jo 

with: 

^5 = 2||C/o||L„ + (2i^ + ^)ll^olli,™ + (4+ ^)||VC/o||L, +4X||V0o||i,™ 
+ /„^(2+^)||F(.)|||„d. 

Ci=^ 2max(l+2si,i+i5Lz^,„iax(l,a2|X7-/3|) + ^(7^7 + /3)) 

With Gronwall lemma we get: 

Jo 
thus 



< eC^*(C2||Xo||t„ + C3|lVXo|li,„ + C4 /o^ |li^(s)|li™ ds) 



with 



D 



4 4r)2 8 

C2 = 2 + -— , C3 = 4 + et C4 = 2 + - 



2.2 Estimation of the non linear terms 

We will now estimate the non linear terms of equation (fl]) . The following propo- 
sition holds true: 

Proposition 3 Let m > 2 be an integer, Xi and X2 elements ofU™^^. Let us 
define: 

Fl = {Ui.V2)u2+WidzU2 
F2 = {Ui.V2)v2 + WidzV2 
F3 = {Ui.V2)O2+Wid,02 

then we have, for alii £ {1,2,3}; 

I!^.||2,™ < C{\\X,\\2^„, + a'\\\/X,\\2,n)\\^X,\\2,,n\\^X2\\l^ 

where C is a constant (independent of a) of order 1. 
If Xi — X2 then we have, for all i G {1,2,3}; 

\mln. < C,{\\X,\\1„,+ \\S/X,\\IJ' 

where C5 = C5(a) is a constant such that C5 — C5 + C'^a^ (where C5 and C5 
are constants of order 1). 
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Proof. First let us address the term uidxU2- 






(16) 



because iJ™ is an algebra for to > 2. 

We estimate now sup^g(o,a) ll'"i('2^)llm- We first write an estimation in space 

dimension 1 for ip{z) G H^{0,a) with (p{0) = 0: 

= j;2^{z')dMz')dz' 

< 2{j;^\^{z)\'dzj;\dMz)\'dzf' 



Similarly we can write an estimation for ||ui(z)||^: 



(17) 



Therefore (16 1 and (17) give us 



\\uidxU2\\i,„ < 2\\ui\\2,m\\d,Ui\\2,mJ^^o\\dxU2{z)\\'^dz 
< 2\\Xi\\2,m\\S/Xi\\2,m\\VX2\\l^ 



And if ui = U2 we have: 



hl5,Ui||2 < 2||Xi||2,™||VXi||3 



< (ll^illi,™+l|VXi|||J 



We establish the same inequalities for the terms uidxU2, Uidx02, vidyUi, vidyV2 
and Vidy02- 



Let us now focus on widzU2. As in (16 1, we have: 

\\widzU2\\i,„< / \\wi{z)\\l^\\dzU2{z)\\l^dz 



We have also, as in (17) 



Thus: 



|wi(z)||^ < 2||w;i||2,,„||92-u;i||2,., 



\hldzU2\\ljn < 2||wi||2,m||92Wi||2,m|!92U2|!i; 



< C\\wi\\2,m\\^X,\\2^ra\\yX2\\l, 



Then we can bound ||wi||2,m as in (12 I: 



kill!,. 



= E|„|<„ /o" \d^y lo d.u,{z') + dyV,{z') dzf dz 

< Eic.\<m lo^dzj;; 2\dSydxU,{z)\^ + 2\d^ydyv,{z)\^ dz 

< a^i\\d.u^\\l^ + \\dyv,\\lj 

< a'l|VXi||2 



therefore we get: 



\\widzU2\ 



2,m 



< Ca^VX^\\lJ\VX2\\l 
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And if u\ = U2 we have: 

\\wid,uA\l^ < CaHVXAlt 



4 
2,m 



< Ca^{\\X 



2 

1||2,; 



IVXil 



We have the same estimates for widzV2 et Widz02, and that concludes the proof 
of proposition l3] 



D 



2.3 End of the proof 

We will now construct a solution of the non linear equation (jTl). Let us first 

introduce some notations. 

Let L be an operator defined as follows: 



L{X,po) = 

And let us define F{Xi,X2)hy 

F{X^,X2)^ 



dfU — Au — av + dxPo + /q /3d.j.,9 dz 
dtv — Av + au + dyPo + Jq l3dy9 dz 
dtO -A0 + -fw 



-{Ui.V2)u2 - WidzU2 
-{Ui.W2)v2 -WidzV2 
-{Ui.\/2)02-Widze2 



We now define N, which is the square of a norm: 



N{X{t)) = \\X{t)\\l„ 



\dtXis)\\l^ds 



We are then given an initial condition Xq £ U"'^^^. Let us define the sequence 
(X",Pq) in the following way: 

X^{t, X, y, z) = Xq{x, y, z), V(i, x, y, z) G M.^ x Q. 

{ i(X"+\pJ,'+i)=F(X",X") pourn>0 

We now verify that the sequence (X",Pq) is well defined for all t and that 
7V(X"(i)) is finite for all n and all t: we denote by N° the value iV(X") = 
ll-'^oll^n.+i- Let us assume that X" is well defined with a finite norm N'^/'^ for 
all t. Thanks to propositions [2] and [sj we get existence and unicity of X"+^ and 
moreover: 

N{X-+\t)) < Coe^^\N^ + f;^\\F{X-,X-)\\l^ds) 

< Coe^^\N' + C,J^\\X-ru^^,ds) 

< Coe^^\N" + C5isup,g[o^,] N{X^{s))^) 

< Cee^^\N° + isup,g[o,,] N{X^{s))^) 

< CXD 

Now, let us choose t* < (4Cge^"-^^* N'^)~^, we then have by recurrence 

sup iV(X"(t)) <2C6e^'''N° 
te[o,f] 
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In fact, we have iVo < 2(766*^^* N'^ (after increasing Cg if necessary) and then 
by recurrence: 

sup,g[o,t-]^^(^"+Hi)) < C6e^i**(iVO + rsup,e[o,,]A^(X»(s))2) 
< C6e^i*'(iV0 + r(2C6e^i*'iV0)2) 

The sequence (X") is therefore bounded for N uniformly in t e [0,f*]. 
To pass to the hniit, we write the equation verified by X"+^ — X": 

L(X"+i - X") = i^(X", X") ~ F(X"-\ X"-i), X"+i - X"|t=o = 

Thanks to propositions [2] and |3] we obtain: 

sup,g[o^,.]7V(X"+i-X") < C/o**||F(X",X")-nX"-i,X"-i)|l|„dt 

= C/o* ||ii^(X"-X"-\X" + X"-i) 

+ iF(X" + X"-i, X" - X"-i)||2_„ dt 

< C/o* 7V(X" - X"-i)(7V(X") + A^(X"-i)) dt 

< CN° J^' N{X" -X"-^)dt 

< CiVOrsup,g[o,,.]iV(X"-X"-i) 

< (CiVOr)"sup,g[o,*l^(^'-^°) 

We can again decrease t* (if necessary) such that CN^t* < 1, consequently 
(X") is a Cauchy sequence with respect to N, and thus it converges to X 
strongly in C{[0,t*],U"'+^); ^jX" converges to dtX strongly in L'^{0,t*; LlH^y) 
and inequality ([3]) is true. 

We can then pass to the limit into the variational formulation of the equation, 
and reintroduce the pressure, thanks to very classical methods (see for example 
[8j and the further paper by Temam and Ziane [13 J. and this concludes the 
proof of theorem [l] 

3 Existence of an optimal control 

In this section we prove theorem [2] 

The proof uses the minimizing sequences method, in two steps: first we prove 

the convergence of the observation term, and then we pass to limit in the state 

equation. 

3.1 Convergence of the observation term 

Let (Xq) be a minimizing sequence for J^: 

J{X^^) ^ inf J{X^) = inf J°(Xo) + ||Xo||^™+, 

Then (XJ) is bounded in ^'"+1 and in (7i™+^)'^, and converges to X^ weakly 
in ^™+i and in (7^™+^)'^. The core of the proof is the convergence of the ob- 
servation term J°{Xl^) = 5l|C"(ii) - ^11^ to \U*{ti) - df where ^ et ^ are 
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the Lagrangian trajectories associated to Xq and Xq. Indeed, we will then get 
the strong convergence of Xq to Xq in W"^^ and in (7^™+^)-^ and we will easily 
verify that Xq is a mininiizer of S . 

Let X" be the solution of equations (ITl and (pi associated with Xq. We 
choose t* according to theorem IT] {t* depends on the norm of Xq and we 
know that the norm of Xq is bounded, thus we can find t* satisfying for 
every X"). This proposition states that the sequence {X"^) is bounded (uni- 
formly in n) in C([0, t*], (Ti"*"^^)^) and that the sequence dtX"^ is uniformly 
bounded in L^(0, t*\ {LlH"^)^). We then prove that {X") is uniformly bounded 
inCV2([0,r],(L27?™)3): 



\\x-ih)-x-{h)h.^ = \\J,'^dtX-h,, 



< 



*^iiaT"i|2 u/2 



vh^^iu:: \\dtx-\\ij 



where C does not depend on n or any U. Thus we get: 

(X") is uniformly bounded in C([0, t*],{n"'+^f) 
(X") is uniformly bounded in C^/^{[0, t*],{LlH^yf) 
=> (X") is uniformly bounded in C'^^^{[0,t*],{[n"^+\ LlH^y]e)^) 

for all 6 e [0,1]. 

We now prove that there exists a small 9 > such that, for 5 > small enough, 

the space [7i™+^, L^i?™]^ is compact in the following space C^: 

C^ = {ue Hl/'^+^Hl+^, periodic in x, y, u = on {z == 0, z == a} x T^} 

For (5 > 0, £'' is a subset of the space of the functions continuous in z and 

Lipschitz in {x,y). 

We now describe the Hilbert interpolation [TC™^^ , LlH"y]g, thanks to Fourier 

series for {x,y) € T^ and sine series in z S (0,a). Let C — {£,-,7]) £ M^ bg the 

Fourier variable associated with the Fourier series on the torus, and k e N* the 

Fourier variable associated with sine series on (0, a). 

We then have: 

^(a;,y,z)eL2i/™ ^ Wi{C,k) ipiCn) G e{Z^ xW) 
ip{x,y,z)en"'+^ ^ W2{C,K)ip{C,K)ef{Z^xW) 

where Wi and W2 are weights defined as follows: 

wi{c,k) = (1 + icn^ 

W^2(C,^) = (l + |C|2)^ + (l + |cn?(l + |«:|2)i 
For all > small enough we get: 

^ £ [H"'+\ LlH^y]g ^ W^Wl-'^ip e 1^{I? X N*) 
with: 

WfWl-' = (l + |CnT^((l+|C|2)=^+(l+|C|2)l(l+|A.|2)l)(l-«) 

= (i + icp)^((i + icp)^ + (i + kn^)(^-''^ 
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For (5 > small enough, we also have: 

where: 

o 2 + 5 o '5 + 1/2 

To establish that the injection [7i™+^, L^iJ™]^ is compact (for 6 > Q small 
enough) in C^ (for (5 > small enough) it is sufficient to find 5 and 9 such that: 

lin, MCi^^^ .0 

|C|,|«H+ooW^l(C,K)«t¥2(C,«)l-' 

We have: 



WfW^-'^ (1 + |CP)^((1 + |CP)5 + (1 + |K|2)^)(i-e) 

(l + |C|2)l(l + |,|2)^±^ 



< 



< 



((i + |CP)5 + (i + l«p)5)(i-«) 
i(i + |CP)f + 1(1 + ^2) ^^^ 

((l + |C|2)5+(l + |K|2)5)(l-e) 



using that m > 2 and ab < - — h — for a, 6 positive and p = 3, q = |. Then we 
get: 

^S /w I *ns 3i5 1-e ,^ , ,„, 35 + 3/2 1-8 

wfi^ < (i+icn— + (i+N^)^-- 

For 6 and (S positive and small enough we have: 

3(5 1-6* 1 3(5 + 3/2 1-61 1 

~' -:^ et ^ — ~ — - 



2 2 2 

Then we can find 9 and S positive and small enough such that the right hand side 
of equation ( |3.1[ ) converges to 0. We have finally shown that (X") is uniformly 
bounded in C^([0,r], ([T^^+i,^^^™]^)^) with [K"'+\ LlH^y]e compact in 
£■'. Then (X") converges strongly to X* in C{[0,t*],{C^f). And finally we 
have: 

Then we prove that (Xq) converges strongly to Xq in U"^'^^: 

J{X^) - ||Xo*||^,„^,+^°(Xo*) 

< lim||^o"l!^™+i+l™^°W) 

< infx„g;^™ + l ^(Xo) 

thus J{X^) - infx„6w..+i^(Xo), then WX^Wu-^-,! ^ |lXo*||;^„.+i, therefore 
(Xq ) converges strongly to Xq in W^~^^, and Xq is a minimizer of J^. 
The strong convergence of (X") in C([0,i*], (C^)^) implies also: 

X*{t^O)^X* 
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3.2 Passage to the limit in the state equation 

In this section we prove that the hniit X* satisfies equation ([T| and boundary 
conditions (pi . To do so, we first pass to the hmit in the weak formulation of the 
equation. Let X' = {u',v',6') € (l?((0,t*) x fi))'^ be a test function satisfying 
conditions (pj). We now pass to the hmit in the five following terms: 

Ti = // dtu^-u' + dtv^-v' + RdtO^-e' 

Ts = j7(C/".V2w" + w"92u")w' + ([/". V2w" + u;"5;,d")u' 
+K{U'\V20" + w"9^6i")6'' 

Ti. As we proved that (9tX") is uniformly bounded in L'^{0,t*; {LlH^yf), it 
converges thus weakly in L^(0,i*; (L^H™)^). We can prove easily that its limit 
is dfX* , e. g., for ^^m": 



dtu^'u' = - / / u"dtu -^11 udtu = / / dtuu' 

Then we get for Ti: 

Ti -^ J J dtu*u' + dtv*v' + Kdte*e' 
T2. A space integration by parts gives: 



To the limit we obtain: 

T2^ - fJiyu*Au' - fJiyv*Av' - K Jfi'e*Ae' 

T3. We easily have: 

T3 -^ JJ —av*u' + au*v' 

T4. The first part is easy. For the second we have: 

= Jj{d,u- + dyv-)[j^,^j'{z')dz'] 

thanks to Fubini's equahty. As (X") converges strongly in i^(0,i*; (L^i/™)'^) 
with 771 > 2, we can pass to the limit, using Fubini again and posing w* — 
J^{d^u*{z')+dyV*iz'))dz': 

Ti^ ~f3jj9*w' -KjJJw*9' 

T5. Let us consider the following part (the other parts can be evaluated simi- 
larly) : 
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We know that dzU^ — dzU* converges weakly to in L^{0, t*; L'IL'^ ). We also 
have w* e L^{0,t*;LlH^y) because X* e L^{0,t*;U"'+^). We now prove that 
w*u' G L'^{0,t*;LlLly) as in the proof of proposition pj 

^ It \\u'\\LlH^J\9zU'\\LiH^^J\w*\\l2H^^^ 

< SUPt [\\u'\\l2H^ \\dzU'\\L2H2 ] \\W*\\1, 2H2 

< oo 

Thus Ts^i converges to 0. 
For T5 2 we have: 



2 



I JJ{w- - w*)dzu-uf < llu;" - ^«11i.((o,*.)xa)ll5^«""'lli^((o,r)xn) 

^ r^llvn \^*||2 IIQ ,,n||2 

^ '-'ll^ ^ 11^2(0 t..i2j^l^)||O^U ||^2(o^f..i2^2^) 

^ /^ll vn Tr*l|2 II V"l|2 

^ '-'ll^ -^ 11^2(0^4. .^2^1^)11-^ llL2(o,t.;W-,^+l) 

SO T5^2 converges also to 0. This gives the expected limit for T5: 

T5 -^ JJ{U*.V2U* + w*dzu*)u' + (C/*.V2W* + w*a^w*)w' 
+ii:(C/*.V26'* + u;*a20*)6l' 

Then the limit X* satisfies the following equation, for all X' e (23>((0, t*) x fi))^ 
satisfying ([2]): 

- //(j/u* Am' + z^w* Au' + Kve*M') 
+ ni-oi'"*^' + otu*v') 

+ n{-pe*w' - if7[/;(a.^i*(^') + dyv*{z')) dz']e') 

- // ( (C/*.V2U* + u;*92U*)u' + {U*.V2V* + w*dzV*)v' 

+K{U*.V20* + w*dz0*)e' ) 

We now prove that this defines a linear functional which is continuous on 
L^(0,i*; (-ffQ(ri))'^. It is obvious that the terms involving linearly X* and its 
derivates define a linear functional. It remains to study terms involving X* non 
linearly. Let us define: 

0(w') — {u*dxU* + w*dzU*)u' 



We will prove that is continuous on i^(0, t*; L'^{fl)) and therefore on £^(0, t*; iJg (fi)). 
As we have X* E C(0, t* ,U™'^^), every space derivative of u* , and also w* and 
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dzW* , are in C{0,t* ,L^H^ ). Then we have, as in the proof of proposition pi 



JJ\u*d^u* +w*dzu*\'^ = J^J\u*dxU* + w*dzU*\\ 



L2 



< 



Cj^J\u*d,u*+w*d,u*\\l, 



< Cj^Ju*\\l,\\d, 



2 



< 



'•'^ xy xy xy xy 

C It\\^*\\L^.mJ\dzU*\\L2H^ \\d^U*\\jj2 

'■ ^ •'■ y ^ -i- y xy 

+ C Jt\\w*\\L2m \\dzW*\\L2H2 \\dzU*\\jj2 



< csupj|x*(i)r 



w™+l 



Thus </) is a continuous hnear functional on L^{0,t*;HQ{il)). Similarly we deal 
with every other non linear term in X* . 



Finally, formula ( 18 ) defines a linear functional, which is continuous on L^ (0, t* ; Hq (fl))^ 
and vanishes on the elements of i^(0, t*; {Hq{U))^ satisfying the boundary con- 
ditions (pi). As in |T5] and in the theory of Navier-Stokes equations (e.g., see 
[71 [TT]), we can reintroduce the pressure: there exists p* G l?'(0,i*, L^(ri)) 
such that {X* ,p*) is a solution of the Primitive Equations (fTl) with boundary 
conditions ^. 

4 Numerical experiments 

In this section, we briefly present some illustrative numerical results. A detailed 
description of the numerical setup and more in-depth results can be found in 

m- 

4.1 Numerical setup 

This problem was addressed with a realistic state-of-the-art Primitive Equations 
ocean model, namely OPA code, developed by LODYC (see [S]). The model is 
set-up in a classical double-gyre wind-driven configuration: it is representative 
of mid-latitude ocean circulation, where a non-linear and non-stationary jet- 
stream (such as the Gulf Stream) develops at the convergence of the subpolar 
gyre and the subtropical gyre. 

The control optimal problem that we solved numerically is very similar to 
the one presented before, the regularization term in the cost function (pi being 
different: 



^(^o) - hJ:-Lij:li\uu)-dir + ^\\xo-x,\\i 



(19) 



The observation term J° is the same, except we have M floats drifting and N 
time-sampling of their positions. The background term J^ involved a so-called 
background state X^,, containing a priori information (such as climatology, or 
results of a previous forecast) and a background error covariance matrix B to 
define the B-norm: 

\\X\W = ^XB-^X 
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observation cost function 




cost function gradient 




Figure 1 : Observation term of the cost function (top) and gradient of the total 
cost function (bottom) as a function of the iterations number. The hnearization 
re-occurs every ten iterations. 
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Figure 2: Kinetic energy of the true state, background (no assimilation) and 
assimilated state at the surface after thirty days; the x-coordinate represent 
longitude grid-points and y-coordinate are latitude grid-points. 
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Figure 3: Relative RMS errors for u (top) and v (bottom) on the whole grid as 
a function of time for the background and the assimilated state, as a function 
of the time. 
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4.2 Four-dimensional variational data assimilation 



The cost function (19 1 is minimized thanks to an iterative process involving its 
gradient. Said gradient is computed thanks to the adjoint of the ocean model, 
and also the adjoint of the observation operator ^. This method, called four- 
dimensional variational data assimilation, has been introduced in meteorology 
by Le Dimet (see [HIS]) and is based on [6]. Although this method has been 
introduced for linear models and observation operators, it can be extended to 
"reasonably" non- linear cases: in [1 the authors proposed an incremental ap- 
proach, which consists in linearizing the operators around a given trajectory 
and then proceed to minimize the quadratic cost function. Once the minimum 
is reached, the operators are once again linearized around it, and the new cost 
function is minimized, and so on. Figure [l] presents the decrease of the cost func- 
tion and its gradient as a function of the iteration number. Every ten iterations, 
the linearized operators and their adjoints are updated. 

4.3 Illustration 

The test-experiment presented here involves M = 1000 floats, drifting at 1000 
meters depth, whose positions are sampled every day during ten days. In order 
to achieve convergence, we had to assume that we have a priori information on 
temperature and salinity, so that the assimilation process aims to reconstruct 
the velocities, from the positions information. We performed identical twin ex- 
periments: a given output of the ocean model is called "true state" and is used 
to generate observations. Then the iterative process is initialized with the back- 
ground (which is equal to the true state except it has wrong velocities), and 
the observations are assimilated. We then get an "assimilated stated" which 
should be close to the true state. The experiment presented here consists of 
three assimilation process on three successive ten-days periods. 
Figure [2] presents the square root of the kinetic energy at time T = 30 days, 
at the surface, for the true state (reference), the background (no assimilation) 
and the assimilated state. The assimilation process reconstructed very well the 
true state. Let us note that the horizontal velocity field presented in Figure [2] is 
the surface one, whereas the floats drift at 1000 meters depth. The assimilation 
thus transfered information to every vertical level and not only the 1000 meters 
deep one. 

Figure Is] presents the evolution of the relative RMS error (with respect to the 
true state) over the 30 days time window. For any velocity fleld {u{x, y, z, t), v{x, y, z, t)), 
such relative RMS error for u is a function of time, computed as follows: 

_ ^ J^\ut{x,y,z,t) -u{x,y,z,t)\'^ dxdydz \in 
^ JQ_\ut{xTy^z^t)\'^ dxdydz I 

where Ut is the true velocity. Similarly we can compute E(v\€). We then can 
compare the errors for the background E(u\y\ t) and the assimilated state £{ua] t). 
We can see that after thirty days the error has been divided by 3 for u and 4 
for V. 
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